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Abstract. We introduce relaxed scheduling as a paradigm for mesh 
maintenance and demonstrate its applicability to triangulating a skin 
surface in R 3 . 
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1 Introduction 

In this paper, we describe a relaxed scheduling paradigm for operations that 
maintain the mesh of a deforming surface. We prove the correctness of this 
paradigm for skin surfaces. 

Background. In 1999, Edelsbrunner [5] showed how a finite collection of spheres 
or weighted points can be used to construct a ^-continuous surface in R 3 . It 
is referred to as the skin or the skin surface of the collection. If the spheres 
represent the atoms of a molecule then the appearance of that surface is similar 
to the molecular surface used in structural biology [2,8]. The two differ in a 
number of details, one being that the former uses hyperboloids to blend between 
sphere patches while the latter uses tori. The skin surface is not C 2 -continuous, 
but its maximum normal curvature, k, is continuous. This property is exploited 
by Cheng et al. [1], who describe an algorithm that constructs a triangular mesh 
representing the skin surface. In this mesh, the sizes of edges and triangles are 
inversely proportional to the maximum normal curvature. The main idea of the 
algorithm is to maintain the mesh while gradually growing the skin surface to 
the desired shape, as illustrated in Figure 1. The algorithm thus reduces the 
construction to a sequence of restructuring operations. There are edge flips, 
which maintain the mesh as the restricted Delaunay triangulation of its vertices, 
edge contractions and vertex insertions, which maintain a sampling whose local 
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Fig. 1. The mesh is maintained as the surface on the left grows into that on the right 



density is proportional to the maximum normal curvature, and metamorphoses, 
which adjust the mesh connectivity to reflect changes in the surface topology. 
Some of these operations are easier to schedule than others, and the most difficult 
ones are the edge contractions and vertex insertions. They depend on how the 
sampled points move with the surface as it deforms. The quality of the mesh is 
guaranteed by maintaining size constraints for all edges and triangles. When an 
edge gets too short we contract it, and when a triangle gets too large we insert 
a point near its circumcenter. Both events can be recognized by finding roots 
of fairly involved functions. Scheduling edge contractions and vertex insertions 
thus becomes a bottleneck, both in terms of the robustness and the running time 
of the algorithm. 

Result. In this paper, we study how fast edges and triangles vary their size, and 
we use that knowledge to schedule these elements in a relaxed fashion, fn other 
words, we do not determine when exactly an element violates its size constraint, 
but we catch it before the violation happens. Of course, the danger is now that 
we either update perfectly well-shaped elements or we waste time by checking 
elements unnecessarily often. To avoid the former, we introduce intervals or gray 
zones in which the shapes of the elements are neither good nor unacceptably bad. 
To avoid unnecessarily frequent checking, we prove lower bounds on how long 
an element stays in the gray zone before its shape becomes unacceptably bad. 
These bounds are different for edges and for triangles. Consider first an edge uv. 
Let R = \\u — v\\/2 be its half-length and g = 1/ max{«;(u), k(v)} the smaller 
radius of curvature at its endpoints. We use judiciously chosen constants C, Q 
and Qi and call the edge 

acceptable "1 ( C/Qo < R/ g, 

borderline I if I CjQ x < R/g < C/Q , 
unacceptable) [ R/g<C/Q\. 

The middle interval is what we called the gray zone above. Assuming uv is 
acceptable, we prove it will not become unacceptable within a time interval of 
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duration At = {29 - 9 2 )g 2 , where 

= RQi - Cg 
RQi + Cg 

In the worst case, R is barely larger than Cg/Qo, so we have 9 > (Qi — Qo)/(Qi + 
Qo) as a worst case bound. We will see that C = 0.06, Qo = 1-6 and Qi = 2.3 
are feasible choices for the constants, and that for these we get 9 > 0.179. . . 
and At I ' g 2 > 0.326 . . .. Consider next a triangle uvw. Let R be the radius of its 
circumcircle, and g = 1/ max{d(ti), n(v), n(w)} the smallest radius of curvature 
at its vertices. We call uvw 

acceptable *| ( R/g<CQo, 

borderline > if < CQ < R/g < CQi, 
unacceptable J [ CQi < Rj g. 

Assuming uvw is acceptable, we prove it will not become unacceptable within a 
time interval of duration At = (29 — 9 2 )g 2 , where 

9 = 1- VR/(CQ ie ). 

In the worst case, R is barely smaller than CQog, so we have 9 > 1 - VQa/Qi. 
For the above values of C, Qo and Qi, this gives 9 > 0.086 . . . and At/ g 2 > 
0.165.... It seems that triangles can get out of shape about twice as fast as 
edges, but we do not know whether this is really the case because our bounds 
are not tight. 

Outline. Section 2 reviews skin surfaces and the dynamic triangulation algo- 
rithm. Section 3 introduces relaxed scheduling as a paradigm to keep track of 
moving or deforming data. Section 4 analyzes the local distortion within the 
mesh and derives the formulas needed for the relaxed scheduling paradigm. Sec- 
tion 5 concludes the paper. 

2 Preliminaries 

In this section, we introduce the necessary background from [5], where skin 
surfaces were originally defined, and from [1], where the meshing algorithm for 
deforming skin surfaces was described. 

Skin surfaces. We write Si = (zi,ri) for the sphere with center Zi € K 3 and 
radius and think of it as the zero-set of the weighted square distance function 
ft : R 3 R defined by f t {x) = \\x - z t \\ 2 - r 2 . The square radius is a real 
number and the radius is either a non-negative real or a non-negative multiple 
of the imaginary unit. We know how to add functions and how to multiply them 
by scalars. For example, if we have a finite collection of spheres Si and scalars 
J2li = 1 then J2"fifi 1S a g am a weighted square distance function, and we 
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denote by S = ^"fiSi the sphere that defines it. The convex hull of the Si is 
the set of such spheres obtained using only non-negative scalars: 

T = {±2^ l]>>* = land ^ °' Vl } • 

We also shrink spheres and write \f~S~ = (z,r/\/2), which is the zero-set of 
2/— f(z). The skin surface defined by the Si is then the envelope of the spheres 
in the convex hull, all scaled down by a factor 1/V2, and we write this as 
F = env \[? "'. Equivalently, it is the zero-set of the point-wise minimum over 
all functions 2/ — f(z), over all S G T , where / is the weighted square dis- 
tance function defined by S. At first glance, this might seem like an unwieldy 
surface, but we can completely describe it as a collection of quadratic patches 
obtained by decomposing the surface with what we call the mixed complex. Its 
cells are Minkowski sums of Voronoi vertices, edges, polygons and polyhedra with 
their dually corresponding Dclaunay tctrahcdra, triangles, edges and vertices, all 
scaled down by a factor 1/2. Instead of formally describing this construction, we 
illustrate it with a two-dimensional example in Figure 2. Depending on the di- 




Fig. 2. The mixed complex decomposes the skin curve and the area it bounds 

mension of the contributing Delaunay simplex, we have four types of mixed cells. 
Because of symmetry, we have only two types of surface patches, namely pieces of 
spheres and of hyperboloids of revolution, which we frequently put in Standard 
Form: 

+ + el = i? 2 , (i) 
g+al-a = ±r 2 , (2) 

where the plus sign gives the one-sheeted hyperboloid and the minus sign gives 
the two-sheeted hyperboloid. 
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Meshing. The meshing algorithm triangulates the skin surface using edges and 
triangles whose sizes adapt to the local curvature. Let us be more specific. At 
any point x € F, let k(x) be the maximum normal curvature at x. In contrast 
to other notions of curvature, k is continuous over the skin surface and thus 
amenable to controlling the local size of the mesh. Call g(x) = 1/k(x) the local 
length scale at x. The vertices of the mesh are points on the surface. For an edge 
uv, let R uv = \\u — v\\/2 be half its length, and for a triangle uvw, let R uvw be 
the radius of its circumcircle. The algorithm obeys the Lower and Upper Size 
Bounds that require edges not be too short and triangles not be too large: 

[L] Ruv/Quv > C/Q for every edge uv, and 
[U] Ruvw/Quvw < CQ for every triangle uvw, 

where g uv is the larger of g(u) and g(v), g uvw is the minimum of g(u), g(v) and 
g(w), and C and Q are judiciously chosen positive constants. 

The particular algorithm we consider in this paper is dynamic, in the sense 
that it maintains the mesh while the surface deforms. We can use this algorithm 
to construct a mesh by starting with the empty surface and growing it into the 
desired shape. This is precisely the scenario in which our results apply. To model 
the growth process, we use a time parameter and let Si{t) — (z i} y/rf + t) be 
the i-th sphere at time t £ R. We start at t = — oo, at which time all radii 
are imaginary and the surface is empty, and we end at t = 0, at which time 
the surface has the desired shape. This particular growth model is amenable to 
efficient computation because it does not affect the mixed complex, which stays 
the same at all times. Each patch of the surface sweeps out its mixed cell. At any 
moment, we have a collection of points sampled on the surface, and the mesh is 
the restricted Delaunay triangulation of these points, as defined in [4,7]. Given 
the surface and the points, this triangulation is unique. As the surface deforms, 
we move the points with it and update the mesh as required. From global and 
less frequent to local and more frequent these operations are: 

1 . topology changes that affect the local and global connectivity of the surface 
and the mesh, 

2. edge contractions and vertex insertions that locally remove or add points to 
coarsen or refine the mesh, and 

3. edge flips that locally adjust the mesh without affecting the point distribu- 
tion or the surface topology. 

For the particular growth model introduced above, the topology changes are 
easily predicted using the filtration of alpha complexes as described in [6]. To 
predict where and when we need to coarsen or refine the mesh is more difficult 
and depends on how the points move to follow the deforming surface. This is 
the topic of this paper and will be discussed in detail in the subsequent sections. 
Finally, edge flips are relatively robust operations, which can be performed in a 
lazy manner, without any sophisticated scheduling mechanism. 

Point motion. To describe the motion of the points sampled on the skin surface, 
it is convenient to consider the trajectory of the surface over time. Note that the 
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i-th sphere at time t is Si (t) = / t ~ (t) . Similarly, the convex combination defined 
by coefficients ji at time t is S(t) = where f = Y1 lift- We can represent 

the skin surface in the same manner by introducing the function g : M 3 — > M 
defined as the point-wise minimum of the functions representing the shrunken 
spheres. More formally g(x) — min{2/(x) — f(z)}, where the minimum is taken 
over all spheres S £ T and z is the center of S. The skin surface at time t is then 
F(t) = <7 _1 (t), so it is appropriate to call the graph of g the trajectory of the 
skin surface. We see that growing the surface in time is equivalent to sweeping 
out its trajectory with a three-dimensional space that moves through time. It is 
natural to let the points sampled on F(t) move normal to the surface. For a point 
x = [£i,£2j S,3] T on a sphere or hyperboloid in Standard Form ^\ + ^±£3 = ±-R 2 , 
the gradient is Vg x — 2[£i, £2, ±£s] T - The point x moves in the direction of the 
gradient with a speed that is inversely proportional to the length. In other words, 
the velocity vector at a point x is 

. _ dx Vg x \7g x 

dt ~ IIV.9,11 2 " 4||x|| 2 ' 

The speed of x is therefore ||±|| = l/(2||x||). The implementation of the relaxed 
scheduling paradigm crucially depends on the properties of this motion. We use 
the remainder of this section to describe a symmetry property of the velocity 
vectors that is instrumental in the analysis of the motion. Consider two mixed 
cells that share a common face. The Standard Forms of the two corresponding 
surface patches differ by a single sign, and so do the gradients. If we reflect 
points in one cell across the plane of the common face into the other cell then we 
preserve the velocity vector, as illustrated in Figure 3. We use this observation 




Fig. 3. Velocity vectors of a shrinking circle on the left and of a hyperbola on the right. 
The right portion of the edge uv is reflected across face shared by the two mixed cells 



about adjacent mixed cells to relate the velocity vectors of points in possibly non- 
adjacent cells. Consider points u and v and let x\, X2, ■ ■ ■ , Xk be the intersection 
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points with faces of mixed cells encountered as we travel along the edge from u 
to v. Starting at i = k, we work backward and reflect the portion of the edge 
beyond x, across the face that contains Xj. In the general case, this portion is a 
polygonal path that leads from x, to the possibly multiply reflected image v of v. 
After k reflections we have a polygonal path from u to the final v. The length of 
the path is equal to the length of the initial edge, and hence \\u — v\\ < \\u — v\\. 
We note that v does not necessarily lie in the mixed cell of u, but its velocity 
vector — which is the same as that of v — is consistent with the family of spheres 
or hyperboloids that sweeps out that mixed cell. In other words, the motion of 
u and v is determined by the same quadratic function. 



3 Relaxed Scheduling 

In this section, we introduce relaxed scheduling as a paradigm for maintaining 
moving or deforming data. It is designed to cope with situations in which the 
precise moment for an update is either not known or too expensive to compute. 

Correctness constraints. In the context of maintaining the triangle mesh of a 
skin surface, we use relaxed scheduling to determine when to contract an edge 
and when to insert a new vertex. Since determining when the size of an edge or 
triangle stops to be acceptable is expensive, we introduce a gray zone between 
acceptability and unacceptability and update an element when we catch it inside 
that gray zone. That this course of action is even conceivable is based on the 
correctness proof of the dynamic skin triangulation algorithm for a range of its 
controlling parameters. The first three conditions defining that range refer to e, 
C and Q. We have seen the latter two before in the formulation of the two Size 
Bounds [L] and [U]: C controls how well the mesh approximates the surface, and 
Q controls the quality of the mesh. Both are related to e, which quantifies the 
sampling density. 

(I) We require < e < sq, where e = 0.279 ... is a root of 2cos(arcsin + 
arcsine) — jz~ — 0. 

(II) Q 2 -4CQ>2. 

It is computationally efficient to select the loosest possible bound for the sam- 
pling density: e = eq. Then we get 5 = 0.166. . . and, as noted in [I], we may 
choose C = 0.08 and Q = 1.65 to satisfy Conditions (I) to (III). Alternatively, 
we may lower C to 0.06 and are then free to pick Q anywhere inside the interval 
from 1.6 to 2.3. The two choices of parameters are marked by a hollow dot and 
a white bar in Figure 5. The last two conditions refer to h, I and m. All three 
parameters control how metamorphoses that add or remove a handle are imple- 
mented. Since the curvature blows up at the point and time of a topology change, 
we use a special and relatively coarse sampling inside spherical neighborhoods of 
such points. Assuming a unit radius of such neighborhoods, we turn the special 
sampling strategy on and off when the skin surface enters and leaves the smaller 
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spherical neighborhood of radius h < 1.0. If the skin enters as a two-sheeted 
hyperboloid we triangulate it using two £-sided pyramids inside the unit sphere 
neighborhood. If it enters as a onc-shcetcd hyperboloid we triangulate it as an 
m-sided drum with a waist. The conditions are stated in terms of the edges ab, 



Fig. 4. The triangulation of a two-sheeted and a one-sheeted hyperboloid inside a unit 
neighborhood sphere around their apices 

be and wx and the triangles abc and vwx, as defined in Figure 4. Their sizes can 
all be expressed in terms of h, I and to, and we refer to [1, Section 10] for the 
formulas. 



Quality buffer. The key technical insight about the dynamic skin triangulation 
algorithm is that we can find constants e, C, h, £, m and Qo < Q\ such that 
Conditions (I) to (V) are satisfied for all Q e [Qo,Qi]. This is illustrated in 
Figure 5, which shows the feasible region of points (C, Q) assuming fixed values 
for e, h, £ and to. Instead of fixing Q and contracting an edge when its size- 
scale ratio reaches C/Q, we suggest to contract the edge any time its ratio is in 
the interval (C/Qi,C/Qo\. After the ratio enters this interval at C/Qq it can 
either leave again at C/Q a or it can get contracted, but it is not allowed to 
reach C/Q\. Vertex insertions are treated symmetrically. Specifically, a triangle 
is removed by adding a vertex near its circumcenter, and this can happen at 
any moment its size-scale ratio is in [CQq, CQi). The ratio can enter and leave 
the interval at CQo, but it is not allowed to reach CQ\. We call (C/Qi, C/Qo] 
and [CQo,CQi) the lower and upper size buffers. The quality of the mesh is 
guaranteed because all edges and triangles satisfy the two Size Bounds [L] and 
[U] for Q — Q\. Symmetrically, the correctness of the triangulation is guaranteed 
because edge contractions and vertex insertions are executed only if the same 
bounds are violated for Q = Qo- 

Early warning. Recall that an edge is borderline iff its size-scale ratio is contained 
in the lower size buffer, and it becomes unacceptable at the moment it reaches 
C/Q\. Similarly, a triangle is borderline iff its size-scale ratio is contained in 




b — c 



(IV) R a b,Rbc,R > C/Q. 

(V) RabcRvwx < mm{(Q,2/Q}Ch. 
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Q (quality) 



Fig. 5. The shaded feasible region of parameter pairs (C,Q) for e = £o, h = 0.993, 
£ — 6 and m = 80. For C = 0.06 this region contains the interval Q £ [1.6, 2.3]. The 
bounding curves are labeled by the corresponding constraints. Redundant constraints 
are not shown 



the upper size buffer, and it becomes unacceptable at the moment it reaches 
CQ\. The relaxed scheduling paradigm depends on an early warning algorithm 
that reports an element before it becomes unacceptable. That algorithm might 
err and produce false positives, but it may not let any element slip by and 
become unacceptable. False positives cost time but do not cause any harm, while 
unacceptable elements compromise the correctness of the meshing algorithm. In 
Figure 6, false positives are marked by hollow dots and deletions are marked by 
filled black dots. All false positive tests of edges are represented by dots above 
the lower size buffer. To get a correct early warning algorithm we just need to 
test each edge often enough so that its size-scale ratio cannot cross the entire 
lower size buffer between two contiguous tests. The symmetric rule applies to 
triangles. Bounds on the amount of time it takes to cross the size buffers will be 
given in Section 4. 

Note that we have selected the parameters to obtain a fairly long interval 
[Qo> Qi]- K is not clear whether or not this is a good idea or whether a shorter 
interval would lead to a more efficient algorithm. An argument for a long interval 
is that the implied large size buffers let us get by with less frequent and therefore 
fewer tests. An argument against a long interval is that large size buffers are 
more likely to cause the deletion of elements that are on their way to better 
health but did not recover fast enough and get caught before they could leave 
the buffers. It might be useful to optimize the length of the intervals through 
experimentations after implementing the relaxed schedule as part of the skin 
triangulation algorithm. 
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Fig. 6. The two buffers are shaded and the two curves are possible developments of size- 
scale ratios for an edge (dashed) and a triangle (dotted). The dots indicate moments 
at which the elements are tested and finally removed 



4 Analysis 

In this section, we derive lower bounds on the amount of time it takes for an edge 
or triangle to change its size by more than some threshold value. From these we 
will derive lower bounds on the time it takes an element to pass through the 
entire size buffer. We begin by studying the motion of a single point. 

Traveling point. We recall that the speed of a point u on the skin surface is 
||u|| = l/(2||w||), assuming we write the patch that contains it in Standard Form. 
The distance traveled by u in a small time interval is therefore maximized if 
it heads straight toward the origin, which for example happens if u lies on a 
shrinking sphere. Starting the motion at point uq, which is the point u at time 
t , we get 

||u|| = yj\\u \\ 2 -(t-t ), (3) 

for the point u at time t. This implies t — to = ||wol| 2 — IM| 2 i so we see that u 
reaches the origin at time t = to + \\uo\\ ■ More generally, we reach the point 
ui = (1 — 9)uq between m and the origin at time ti — t + ||uoll — = 
to + {29 — 2 )\\uo\\ . Since the above analysis assumes the fastest way u can 
possibly travel, this implies that within an interval of duration At = ti — to, the 
point uo cannot travel further than a distance 9g{uo). We use 9 as a convenient 
intermediate quantity that gives us indirect access to the important quantity, 
which is At. 

Recall from the Curvature Variation Lemma of [1] that the difference in 
length scale between two points is at most the Euclidean distance. If that distance 
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is ||uo — mi| < 6q(uq) then the length scale at u\ is between 1—8 and 1 + 9 
times the length scale at uq. It follows that if we travel for a duration At = 
(29 — 9 2 )g 2 (uo), we can change the length scale only by a factor 

i_ fl <££4<i + *. (4) 

The lower bound is tight, and the upper bound cannot be reached because the 
distance 9g 2 (uo) from uo can only be achieved if the length scale shrinks. We 
will also be interested in the integral of l/(2||u|| ), which is again maximized if 
u moves straight toward the origin: 

d* 



J t 2\\u\\ Jt " 



2\\u \\ - (2t-2i ) 

1 



i !m ll M o|| 2 - (h - t ) 
(_ )ln_ 



• I 2 

\ u o\\ 



lnP4 



Mil 



Denoting the above integral by X and choosing t\ — to = {26 — 9 2 )\\uo\\ 2 , as 
before, we have 

x < JW _ g(up) 1 . . 

" ||mi|| q(u x ) -1-9- W 

Edge length variation. Consider two points u and v on the skin surface during 
a time interval [io,ii]- We assume that both points follow their trajectories 
undisturbed by any mesh maintenance operations. Let uo and m be the point 
u at times to and t\ and, similarly, let vo and v\ be the point v at these two 
moments. We prove that if the time interval is short relative to the length scale 
of the points then the distance between them cannot shrink or grow by much. 

Length Lemma. Let g = mm{g(u ), g(v )} and At = h - t = (29 - 9 2 )gl, 
for some < 6 < 1. Then 

1-9 < h-* 1 !! < 1 



||uo-«o|| 1-0 

Proof. The derivative of the distance between points u and v with respect to 
time is 

d||w — u || d||u — v\\du d\\u — v\\dv 
dt du dt dv dt 

(u — v) T 



u — v\ 



(u-v). (6) 



For example if u and v lie on a common sphere patch then g = g(u) = g(v), 
u = ±u/(2g 2 ) and v = ±v/(2g 2 ), which implies 

d||zt — v\\ _^(u — v) T (u — v) ||u — v\\ 

dt ~~ \\u-v\\ 2g 2 ~~ 2g 2 
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We prove below that in the general case, the distance derivative stays between 
these two extremes: 

< d\\u-v\\ \\u-v\\ 
2g 2 ~ dt ~ 2g 2 ' [ ' 

where g — mm{g(u), g(v)}. To get the final result from (6), we divide by \\u — v\\, 
multiply by dt, and use dlnx^dx/xto get 



dt , ... dt 

v - d(ln||u - y|l) ^ v 



Next we integrate over [to,ti] and exponentiate to eliminate the natural loga- 
rithm: 

IK - «o|| 

The claimed pair of inequalities follows from (5) and the observation that the 
upper bound for X cannot be realized when the distance derivative is positive. 
To prove (7) for general points u and v, it suffices to show that the length oiu — v 
is at most ||u — v\\/(2g 2 ). We have seen that this is true if u and v belong to a 
common sphere patch. It is also true if u and v belong to a common hyperboloid 
patch because 

v! v' \\u — v\\ 

2g 2 (u) 2g 2 (v) n ~ 2g 2 ' 

where the primed and unprimed vectors are the same, except that they have a 
different sign in the third coordinate. We need a slightly more elaborate argument 
if u and v do not belong to the same mixed cell. We then reflect v across the 
faces of mixed cells that intersect the edge uv. As described in Section 2, such a 
sequence of reflections does not affect the velocity vector. The distance between 
u and the image v of v under the composition of reflections is at most that 
between u and v. Hence, 

||«-*|| = ||«-*|| < ll^ll, 

as required. ED 

The lower bound in the Length Lemma is tight and realized by points u and 
v on a common sphere patch. 

Shrinking edge. Consider an edge uv, whose half-length at time t is Rq. As 
before, let u and be the points u and v at time t . Let go = min{g(u ), g(vo)}- 
We follow the two points during the time interval [to,ti], whose duration is 
At = t\ — to = (29 — 9 2 )go- The Length Lemma implies that at time t\, the 
length of the edge satisfies 
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Our goal is to choose 9 such that the edge at time t\ is guaranteed to satisfy the 
Lower Size Bound for Q = Q\. Using Ri > (1 — 8)R from (8) and g\ < (1 + 0) go 
from (4), we note that R 1 /g 1 > C/Qi is implied by (l-6)R /(l + 6) > Cg /Qi- 
In other words, 

= RqQi - Cg 

RoQi +Cg K ' 

is sufficiently small. The corresponding time interval during which we can be 
sure that the edge uv does not become unacceptably short has duration At = 
(29 — 9 2 )g 2 . To get a better feeling for what these results mean, let us write 
the half-length of u n v as a multiple of the lower bound in [L] for Q = Qq: 
R = ACgo/Qo with A > 1.0. We then get 9 = (AQ 1 - Q )/(AQ 1 + Q ) and 
At from 9 as before. Table 1 shows the values of 9 and At for a few values of A. 



Table 1. For edges, the values of 9 and At for Qo = 1.6, Qi = 2.3 and a few typical 
values of A 



A 


9 


Atjei 


1.0 


0.179... 


0.326... 


1.5 


0.366 . . . 


0.598... 


2.0 


0.483 . . . 


0.733... 


2.5 


0.564... 


0.810... 


3.0 


0.623 . . . 


0.858... 


3.5 


0.668 . . . 


0.890... 


4.0 


0.703 . . . 


0.912... 



Height variation. Consider a triangle uvw during a time interval [to,ti]- We 
assume that all three points follow their trajectories undisturbed by any mesh 
maintenance operations. Each vertex has a distance to the line spanned by the 
other two vertices, and the height H of uvw is the smallest of the three dis- 
tances. If uv is the longest edge then H = \\w — where w' is the orthogonal 
projection of w onto uv. We prove if the time interval is short relative to the 
length scale at the points then the height cannot shrink or grow by much. To 
state the claim we use indices and 1 for points and heights at times to and t\. 

Height Lemma. Let g = mm{g(u a ), g(v a ), g(w )} and At = ii - to = (29 - 
9 2 )qI, for some < 9 < 1. Then 



Proof. We prove that (7) is also true if we substitute the height H for the length 
of the edge uv: 
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where g = min{ g(u), g(v), g(w)}. The claimed pair of inequalities follows as 
explained in the proof of the Length Lemma. To see (10) note first that the 
height of the triangle is always determined by a vertex and a point on the 
opposite edge, eg. H = ||w — w'\\. Let w' = (1 — X)u + Xv. If u and v belong to 
the same mixed cell then Vg w ' = (1 — A) X7g u + XV g v because the gradient varies 
linearly. Along a moving line segment uv the velocity vectors vary linearly, hence 
w' = (1 — X)ii + Xv. Since the gradients and the velocity vectors at u and v point 
in the same directions, they do the same at w'. The length of the velocity vector 
at w' is at most that of the longer velocity vector at u and v. If w belongs to 
the same mixed cell as w' , this implies 



... \\w-w'\\ _ H 



\w — w \\ < 



2g 2 2g 



2> 



from which (10) follows. If u, v and w do not belong to the same mixed cell 
then we perform reflections, as in the proof of the Length Lemma, and get (10) 
because reflections do not affect velocity vectors. ED 

In the following, we only need the lower bound in the Height Lemma, which 
is tight and is realized points u, v and w on a common sphere patch. 

Expanding triangle. We use both the Length L md the Height Lemma to 

derive a lower bound on the length of time during which a triangle that initially 
satisfied the Upper Size Bound [U] for Q = Q is guaranteed to satisfy the 
same for Q — Q\. We begin by establishing a relation between the circumradius 
R = R U vw of a triangle uvw and its height and edge lengths. Referring to Figure 
7, we let z denote the center of the circumcircle. Assuming uv is the longest of 
the three edges, the height is H = \\w — w'\\ and v and z lie on the same side of 
the line passing through u and w. Let z' be the midpoint of uw and note that the 
angle at z is twice that at v. Zuzw — 2Zz'zw = 2Zuvw. This implies that the 
triangles ww'v and wz'z are similar, and therefore \\z' — w\\/R = H/\\v — w||. It 
follows that the circumradius of uvw is 

\\u — w\\ \\v — w\\ 
R = — 



2H 

There are three ways to write twice the area as the product of an edge length 
and the distance of the third vertex from the line of that edge: ||u — H = 
\\u — w\\ \\v — v'\\ = \\v — w\\ \\u — u'\\. Hence, the circumradius is also 

\\u — v\\ \\u — w\\ \\u — v\\ \\v — w\\ 

R = 



2||u-«'|| 2\\v-v , \\ 

For the remainder of this section, we use indices and 1 for points, heights and 
radii at times to an d t\. The above equations for the circumradius imply 

Ri \\ui - wi|| \\vi - wi|| \\w - w' a \\ 
R \\u - wo\\ \\vo - too || \\wi - will ' 
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Fig. 7. The triangle uvw is similar to wz'z, which implies a relation between the height 
H and the circumradius R 



Assuming Ho = \\wo — w' \\ is the height at time to, we have Hi < \\wi — w[ || at 
time t\. We can therefore use the Length Lemma to bound the first two ratios 
and the Height Lemma to bound the third to get 

We now choose 9 such that a triangle that satisfies [U] for Q = Q at time to 
is guaranteed to satisfy [U] for Q — Qi at time t\. Using i?i < Rq/{1 — 9) 3 
from (11) and (1 — 9)g a < Q\ from (4), we note that R\/ Q\ < CQ\ is implied by 
i?o/(l — 9) 4 < CQiQq. In other words, 

9 = 1 - i/Ro/iCQigo) (12) 

is sufficiently small. It is convenient to write the circumradius of the triangle 
UqVoWo as a fraction of the upper bound in [U]: i?o = CQoQo/A with A > 1.0. 
Then, 9 = 1— \/Qo/(AQi). Table 2 shows the values of 9 and At for a few 
values of A. 



Table 2. For triangles, the values of 9 and At for Qo = 1.6, Qi = 2.3 and a few typical 
values of A 



A 


9 


At/gt 


1.0 


0.086 . . . 


0.165... 


1.5 


0.174... 


0.319... 


2.0 


0.232 . . . 


0.410... 


2.5 


0.273 . . . 


0.472 . . . 


3.0 


0.306 . . . 


0.518... 


3.5 


0.332 . . . 


0.554... 


4.0 


0.354... 


0.583... 
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5 Discussion 

The main contribution of this paper is the introduction of relaxed scheduling 
as a paradigm for maintaining moving or deforming data, and the demonstra- 
tions of its applicability to scheduling edge contractions and vertex insertions 
maintaining skin surfaces. 

Algorithm design. We view the dynamic skin triangulation algorithm, of which 
relaxed scheduling is now a part, as an interesting exercise in rational algorithm 
design. What are the limits for proving meshing algorithms correct? This design 
exercise gives us a glimpse on how complicated meshing problems can be. Per- 
haps more importantly, it illustrates what it might take to prove other meshing 
algorithms correct. We especially highlight the role of constant parameters in the 
algorithm and how they control the algorithm as well as the constructed mesh. 
In our example, the important parameters are C, which controls how closely 
the mesh approximates the surface, and Q, which controls the quality of the 
mesh. The effort of proving the various pieces of the algorithm correct has lead 
to inequalities for these parameters. In other words, we have identified a feasi- 
ble region which is necessary for our proofs and sufficient for the correctness of 
the algorithm. The detailed knowledge of this feasible region has inspired the 
idea of relaxed scheduling, and it was necessary to formulate it in detail and to 
prove its correctness. Many meshing algorithms are based on parameters that 
are fine-tuned in the experimental phase of software design. We suggest that in 
the absence of detailed knowledge of limits, fine-tuning is a necessary activity 
that gropes for a place in the feasible region where correctness is implied. Of 
course, it might happen that this region is empty, but this is usually difficult to 
determine. 

Future work. It is not our intention to criticize work in mesh generation for the 
lack of correctness proofs. Indeed, it would be more appropriate to criticize our 
own work for the lack of generality. Although we laid out a complete algorithm 
for maintaining the mesh of a deforming surface, we are a far cry from being able 
to prove its correctness for any surface other than the skin surface introduced in 
[5] . We have also not been able to extend the algorithm beyond the deformations 
implied by growing the spheres that define the surface. For example, it would be 
desirable to maintain the mesh for deformations used for morphing as described 
in [3]. Generalizing the algorithm to include this application and proving it 
correct may be within reach. 

Another worthwhile task is the implementation of relaxed scheduling as part 
of the dynamic skin algorithm. Are our lower bounds for the necessary At suffi- 
cient to eliminate edge contractions and vertex insertions as a bottleneck of the 
algorithm? Can these lower bounds be improved in any significant manner? Can 
we improve the performance by fine-tuning the parameters, in particular Qq and 
Qi, while staying within the proved feasible region? 



Relaxed Scheduling in Dynamic Skin Triangulation 



17 



Acknowledgments 

The authors thank Robert Bryant for helpful discussions concerning the proof 
of the Length Lemma. 

References 

1. H.-L. Cheng, T. K. Dey, H. Edelsbrunner and J. Sullivan. Dynamic 
skin triangulation. Discrete Comput. Geom. 25 (2001), 525-568. 

2. M. L. Connolly. Analytic molecular surface calculation. J. Appl. Crystallogr. 
6 (1983), 548-558. 

3. H.-L. Cheng, H. Edelsbrunner and P. Fu. Shape space from deformation. 
Comput. Geom. Theory Appl. 19 (2001), 191-204. 

4. L. P. Chew. Guaranteed-quality mesh generation for curved surfaces. In Proc. 
9th Ann. Sympos. Comput. Geom., 1993, 274-280. 

5. H. Edelsbrunner. Deformable smooth surface design. Discrete Comput. 
Geom. 21 (1999), 87-115. 

6. H. Edelsbrunner and E. P. Mucke. Three-dimensional alpha shapes. ACM 
Trans. Graphics 13 (1994), 43-72. 

7. H. Edelsbrunner and N. R. Shah. Triangulating topological spaces. Inter- 
nat. J. Comput. Geom. Appl. 7 (1997), 365-378. 

8. B. Lee and F. M. Richards. The interpretation of protein structures: esti- 
mation of static accessibility. J. Mol. Biol. 55 (1971), 379-400. 

Appendix 

Table 3. Notation for important geometric concepts, functions, variables, and con- 
stants 

fi : R 3 — ► R weighted (square) distance function 

Si = (zi, ri) zero-set of ft; sphere with center z; and radius r» 

T convex hull of spheres Si 

F — env \f~T skin surface 

K ) 6 — V K maximum curvature, length scale 
Qo < Q < Qi constant controlling quality 
e, C, h, £, m additional constants 

g : R 3 -> R point-wise min of the 2/ - f(z) 
F(t) — <? _1 (t) skin surface at time t 



t, 6 time parameter, relative travel distance 

[to,ii] time interval 

At = ti — to duration 

u,u',u point, projection, reflection 

Vg u ,u gradient, velocity vector 



uv,uvw,H,R edge, triangle, height, radius 



